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A generalization of the Cooper pairing mechanism is proposed which allows for a triplet state of 
lower energy. This is achieved by incorporating spin into the canonical commutation relations and 
by modifying the S potential contact interaction. The gap equation contain as solutions both singlet 
and triplet states. It is shown that the triplet state is lower in energy than the singlet state which 
may explain the spin-triplet superconductivity observed in heavy fermion compound UPt3 and in 
Sr2Ru04. 
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Superconductivity [1, 2] is a fascinating subject which 
has been studied vigorously both theoretically as wc^ll 
as experimentally in the past several decades [3, 4]. In 
conventional superconductivity, the electron-phonon in- 
teraction is known to lead to a condensed state of s-wave 
Cooper pairs and this has been experimentally verified. 
However, in certain materials such as the heavy fermion 
compound UPt3 [5, 6] as well as Sr2Ru04 [7, 9], it has 
been observed that the Cooper pairs occur in a triplet 
state. This unconventional superconductivity [10, 11] 
would be compatible with some kind of ferromagnetic 
effect due to the spin in a manner analogous to the spin 
triplet state in superfluid He'^. However, there is no obvi- 
ous way to accommodate this within the framework of the 
conventional description of superconductivity. Any gen- 
eralization of the Cooper pair mechanism should not only 
allow for the existence of both singlet as well as triplet 
states, but should also predict that the triplet state is 
energetically favorable at low energies. 

In the absence of a mechanism to describe the uncon- 
ventional superconductivity within the usual framework 
of quantum mechanics, in this letter we propose a mini- 
mal modification of quantum mechanics that can achieve 
this which can be thought of as a starting point for a 
macroscopic theory for such systems. Basically we mod- 
ify the canonical commutators in such a way that the 
space coordinates and spin become coupled through 

while all the other commutators involving spin and mo- 
mentum remain unchanged. Here is a very small real 
constant parameter with dimensions of length. In fact, 
such a modified algebra was already proposed in [12] 
and studied for the case where the parameter 9 was real. 
Equation (1) merely corresponds to the case where this 
parameter is purely imaginary (namely, letting 9 ^ iO in 
[12]) and in this case the operator Xi ceases to be Her- 
mitian. As a result, the full algebra satisfying Jacobi 



identity also involves the additional relations 

[x^ ^Xj ] — i9 ^ijk'^k^ ^Pj\ — ^^ij 1 

\Xi^ J Sj'\ — ^^ijk^k^ \'^'}} — ^ijk^k' (2) 

It is clear from (1) and (2) that the coordinate opera- 
tors Xi are now matrices in the spin space. Furthermore, 
the complete algebra can also be realized in terms of the 
conventional operators through the change of basis 

Xi = Xi -\- iOSi^ Xi^ = Xi i9Si^ pi = Pi. (3) 

where {xi,Pi,Si) satisfy the standard commutation rela- 
tions. In order to generalize the Cooper pair mechanism, 
we note that given a real (Hermitian) potential V{x), 
a naive generalization to a standard Hamiltonian of the 
form (we set m= 1) 

F(x,p) = ^ + y(x), (4) 

would no longer be Hermitian nor would it be VT sym- 
metric [13] where the transformations of Xi under V and 
T can be seen from (3) to correspond to 

^ ^ t - t - ^ r ^ ^ t - t 

. ^ ■ ' T* ■ ' ^ • ' H' • * T* ■ • ' * - ' 1^1 

U^i 7 Uji , Jji 7 O/^, 7 Jji^ U^i 7 . yU J 

We note here that even a VT symmetric potential V{ix.) 
would not generalize to PT symmetric system and we 
would restrict ourselves to real potential V{x) in this 
work since our goal is to generalize the Cooper pair mech- 
anism. In order to ensure the reality of the spectrum of 
the Hamiltonian, we note that given a potential V{x.), we 
can construct two Hermitian Hamiltonians in a straight- 
forward manner as [14] 

fl-+(x,p) = ^ + ^(y(x) + Ft(x)), (6) 

H_{±,p) = ^+'-{Vi±)-VH^)). (7) 
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Furthermore, we note that if V(— x) = V{ic) then _ff_|- is 
also VT symmetric while if V{—x) = —V{x) then H- 
is VT symmetric as well. A straightforward calculation 
show that the Hamiltonians (6) and (7) can be expressed 
in terms of the standard operators (3) as 



H-{x,p,s) 



Y + l (^(^ + ^^s) + V{x - i9s)) 
^+cos((?s. V)y(x), (8) 
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Y + ^ {v{^ + «^s) - y(x - ies)) 



sin(6's- V)y(x), 



(9) 



where we have used the reality of the potential V{x). We 
note from (8) and (9) that when 0^0 



+ y(x), H_ 
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(10) 



Therefore, we would concentrate on the Hamiltonian 
which would reduce in the limit ^ — >■ to the conventional 
interacting Hamiltonian of the system. Furthermore, in 
spite of the nontrivial mixing of the space and spin vari- 
able, H+ in (8) commutes with s so that the eigenstates 
can be written as product states. 

The generalization of this construction to a many par- 
ticle system can be carried out in a straightforward man- 
ner. For example, let us consider a two particle system 
(which is relevant for our analysis) where the interaction 
depends on the relative coordinate r = xi — X2. First, 
we note that the only combination of the coordinates Xi 
and X2 that commutes with the total spin S = Si + S2 of 
the system is Xi — X2^ and its conjugate. Therefore, the 
Hermitian potential which would also commute with the 
total spin can be written in the form 

y2(ii,x2) = ^ (f(xi - 4) + vHsc, - xt)) 

^^(V{r + teS) + V{r-ieS)), (11) 



where we are assuming that the couplings are real. In 
particular, if we consider a contact potential which is 
responsible for the instability of Fermi sea against small 
attractive interactions between electrons [1, 2, 11], this 
would generalize in our case to a problem of two electrons 
interacting through the "contact" potential 

V2{r, S) = (Sg{r + teS) + <5,(r - z^S)) (12) 

where we have defined a modified delta function as 



|k|= 



|k|=0 



(27r)3 



(13) 



which reduces to the standard delta function in the limit 
^ ^ 0. As a result, we can write the "contact" potential 
as 



|k| = i 



V2iv,S) 



-^1 



|k|=0 



(27r)3 



e'''"'cosh(6lk- S), (14) 



where 7 is a positive constant. We note that the minimal 
modification of the delta potential preserves the qualita- 
tive nature of the electron-phonon interaction (namely, 
the two electrons interact through a local deformation of 
the lattice), but introduces a finite range of interaction 
of the order 20. 

Since the potential depends on the relative coordinate, 
the energy levels are determined from the Schrodinger 
equation for the reduced system 

(-V^ + V2{v, S))*(*'™)(r) = E(«.'")*(«.™)(r), (15) 

where we use the compact notation ^'(«''")(r) = 
^'(r)|s,m) with |s,m) labeling the representation of S. 
In the momentum space 

^-(«.'")(r) = J (fke^^ ''^'^''"'\k), (16) 
we can write (15) as (for |k — k'| < g) 



(k) 



(k2 - 

7 

(k2 - £;(s,m)^ 



j dPk' cosh(6i(k - k') • S)i>('*''")(k') 

j (fk' (cosh(6ik • S) cosh(6ik' • S) + sinh(6'k • S) sinh(6'k' • S)) (k'). (17) 



From (16) we note that the antisymmetry of the two elec- mines (P12 is the exchange operator) 
tron wavefunction under the exchange of particles deter- 

Pi2$(*'"^(k) = (-l)*+i$('^'™)(k), 

(-k) = (-!)''$("''") (k). (18) 



3 



Using these relations, it follows from (17) for the singlet 
state (s = 0) that 



7 



d^k (cosh2(6»k-S))(°'0) 



(27r)3 (k2 - E(ofi)) 
while the triplet state (s = 1) leads to 



= 1, 



(19) 
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J (2^ 



d^k_ (cosh^(6>k-S)-l)(i'' 
)3 (k2 - i;(i.™)) 



(20) 



which can be thought of as the gap equations for the 

singlet and the triplet states respectively. Here (•)('''™) 
represents the expectation value in the state \s,m). We 
can simplify these equations by noting that 

(k • S f = fc2 (1 + (k . 0-i)(k • 0-2)fc-2)/2 = k^Si2, (21) 

where we have identified k = \k\,Si = <Tj/2,i = 1,2 and 
we note that = S12 which leads to 



^^"512, cosh2(6'k-S) = l+sinh2(6'A;)S'i2. (22) 



(k-S)2" 

With a little bit of algebra, we can show that 



dn 



(5l2)(°'°)=0, 



4^ (Si2)(^''") = = ^(1 + ^(1 - (-1)1-1)). (23) 



It is worth noting here that in the triplet state the ex- 
pectation value in the state m = is lower than that in 
the states m = ±1 state which are degenerate and reflect 
a sort of ferromagnetic degeneracy. 

To investigate the Cooper pair mechanism, we convert 
the momentum integrations in (19) and (20) to c;nergy in- 
tegrals through the identification 2e = k^. Furthermore, 
as is well known the electrons which contribute to su- 
perconductivity lie in a thin energy shell near the Fermi 
surface, namely, ep < £ < £f + £c where cp is the Fermi 
energy and Cc <^ ep so that we can restrict the energy 
integration to this small interval. Using (22) and (23) in 
(19), we obtain 



7iV(0) 



de 



(2e-£(o,o)) 



= 1, 



(24) 



where A'"(0) = ^ ^ is the density of states per unit en- 
ergy interval which is assumed to be constant within this 
small range of the integral. This (gap) equation for the 
singlet state coincides with the standard commutative 
case. Therefore, there is no noncommutative correction 
to the energy of the singlet state which has the value 



£1(0,0) =2(eF-ece~'^^ 



(25) 



On the other hand, using (22) and (23) in (20) we 
obtain 



7/9('")iV(0) 



de 



(2e_ £;(i,m)) 



eF+Ec 



27p('")7V(0)6'2 



dee 



(2e-£;(i.™)) 



= 1. 



(26) 



The integration can be carried out in a straightforward 
manner and leads to 



2 / 



1 



In 1 



2e. 



2eF - Ei^'"^) 



£;(!,"») V7iV(0)p('»)6'2 

(27) 

To show that the triplet state has a lower energy than 
the singlet state, let us define the dimensionless ratio 



X = 



£;(o,o) ' 



(28) 



and rewrite (27) as 

C=l + Bxln{l + 



1 



B{1 -x) + e-^ 
where we have defined three dimensionless constants 
2 „ ep . „ A 



(29) 



7^(0)' ec ee ^ 



\ C 



2p(™)6l2ec' 

(30) 

Equation (29) can be solved graphically by looking for 
the simultaneous solutions of the equations 



y = C, y = l + Bx\n[ 1 + 



B{l-x)+e-^ 



, (31) 
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FIG. 1. Simultaneous solution of (31). 

which are plotted in Fig. 1 for different values of the con- 
stants B, C . Wc note that in the weak coupling approx- 
imation A ^ 6.7 and for conventional superconductors 
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B > 10 (recall that ec <C ep) and since is a small pa- 
rameter, the constant C can be chosen to be of the order 
of 10^. The parameters in the plot are chosen keeping 
this in mind. Let us discuss some of the general features 
of the simultaneous solution. First as 9 is taken to be 
smaller and smaller and, therefore, C larger, the solution 
a; — ^ 1 implying that in the limit 6^0, the singlet and 
the triplet solutions become degenerate. However, for 
any finite and small 6, the solution exists for x < I im- 
plying that the triplet solution is lower in energy than the 
singlet solution. Furthermore, from the fact that p^™^ is 
larger for m = ±1 than for m = (see (23)), the constant 
C correspondingly is smaller for m = ±1 than for to = 0. 
This implies that the triplet state with to, = ±1 has the 
lowest (degenerate) energy. The graphical solutions can 
also be understood qualitatively as follows. Since e""^ is 
a very small quantity, when B{1 — a;) < 1, then 

(29) can be written approximately as 



C.^, or ..l-i<l. (32) 



For 6' — > 0, C oo so that a; 1 as we have already 
mentioned. For to = ±1, the value of C is smaller than 
that for TO = and, correspondingly, the value of x would 
be lower for to = ±1. 

The result of this analysis is quite interesting because 
it describes a theoretical mechanism for a triplet pair- 
ing state which is lower in energy than the singlet state 
which may explain the behavior of some of the unconven- 
tional superconductors. Although the precise mechanism 
of why some materials prefer a singlet state rather than 
a triplet state as ground state or vice versa is not under- 
stood yet, our analysis suggests that such an understand- 
ing necessarily requires us to go beyond the conventional 
quantum theory. 

Acknowledgment Two of the authors (JG and FT) 
would like to thank the kind hospitality at the Univer- 
sity of Rochester where part of this work was done. This 
work was supported in part by US DOE Grant num- 
ber DE-FG 02-91ER40685, by MECESUP FSM 0605 and 
by FONDECYT-Chile grant-1095106, 1100777 and Dicyt 
(USACH). 



[1] L. N. Cooper, Phys. Rev. 104, 1189 (1956). 

[2] J. Bardeen, L. N. Cooper and J. R. Schrieffer, Phys. Rev. 

106, 162 (1957). 
[3] P. G. De Gennes, Superconductivity of Metals and Alloys, 

Benjamin, New York (1966) (reprinted, Perseus Books, 

Reading, MA (1999)). 
[4] M. W. Tinkham, Introduction to Superconductivity, 2nd 

edition, McGraw-Hill, New York (1996). 
[5] D. J. Bishop, C. M. Varma, B. Batlogg, E. Bucher, Z. 

Fisk, and J. L. Smith, Phys. Rev. Lett. 53, 1009 (1984) 
[6] H. Ton, Y. Kitaoka, K. Asayama, N. Kinmra, Y. Onuki, 

E. Yamamoto and K. Maezawa, Phys. Rev. Lett. 77, 

1374-1377 (1996) 
[7] T. M. Rice and M. Sigrist, J. Phys. Condens. Matter 7, 

L643 (1995). 

[8] I. Mazin and D. J. Singh, Phys. Rev. Lett. 79, 733 (1997) 



[9] K. Ishida, et al, Nature 396, 658 (1998). 

[10] K. H. Bennemann and J. B. Ketterson in Superconduc- 
tivity; conventional and unconventional superconductors, 
pag. 3, Springer 2008, Vol. I. Edited by K. H. Bennemann 
and J. B. Ketterson. 

[11] See e.g. L. Pitaevskii in Superconducttvity; conventional 
and unconventional superconductors, pag. 27, Springer 
2008, Vol. I. Edited by K. H. Bennemann and J. B. Ket- 
terson. 

[12] H. Falomir, J. Gamboa, J. Lopez-Sarrion, F. Mendez 

and P. A. G. Pisani, Phys. Lett. B 680, 384 (2009) 

[arXiv:0905.0157 [hep-th]]. 
[13] C. M. Bender, S. Boettcher and P. Meisinger, J. Math. 

Phys. 40, 2201 (1999) [arXiv:quant-ph/9809072]. 
[14] D. Bazeia, A. Das, L. Greenwood and L. Losano, Phys. 

Lett. B 673, 283 (2009). 



